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quantum dot Josephson junction 
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We measure the non-dissipative supercurrent in a single InAs self-assembled quantum dot (QD) 
coupled to superconducting leads. The QD occupation is both tuned by a back-gate electrode and 
lateral side-gate. The geometry of the side-gate allows tuning of the QD-lead tunnel coupling in a 
region of constant electron number with appropriate orbital state. Using the side-gate effect we study 
the competition between Kondo correlations and superconducting pairing on the QD, observing a 
decrease in the supercurrent when the Kondo temperature is reduced below the superconducting 
energy gap in qualitative agreement with theoretical predictions. 

PACS numbers: 73.63.Kv, 73.23.Hk, 74.45,+c, 74.50.+r, 85.25.Cp 



Devices which combine the gate tuneability of semi- 
conductor quantum dots (QDs) with non-dissipative su- 
perconducting transport are desirable for very sensitive 
and controllable coherent switching devices 0,0 and the 
study of interplay between Kondo physics and the su- 
perconducting proximity effect 0, 0, 0, HI- The im- 
portant energy scales for the interaction between Kondo 
singlet state and superconductivity are captured in the 
scaling parameter tx = fcs?/f/A determined by the su- 
perconducting gap (A) and Kondo temperature (Tx). 
For 3> 1 the local magnetic moment of the unpaired 
electron spin on the QD is screened by the Kondo ef- 
fect and the ground state of the system is a Kondo sin- 
glet state. In this regime an enhanced supercurrent due 
to the Kondo effect has been predicted theoretically 0- 
However, for tx <C 1 the Kondo state is suppressed by 
the lack of low energy excitations in the superconducting 
energy gap and the system ground state is a degenerate 
(so called magnetic) doublet state. To date experimental 
efforts to elucidate this phase transition have been unable 
to systematically control the physical parameters which 
determine tx due to the limited tuneability of the de- 
vices studied. In the present paper, we demonstrate that 
tx may be smoothly controlled by a side-gate, which is 
placed laterally to an InAs self- assembled QD. The side- 
gate performance is very effective for uncapped InAs self- 
assembled QDs as the lateral confinement tuned by the 
sidegate is weak relative to that in nanowire or nanotube 
devices. When tx is tuned through tx ~ 1 we observe 
a dramatic change in the superconducting transport in- 
dicating the phase transition between Kondo singlet and 
magnetic doublet states. 

To date works on interplay of Kondo and proxim- 
ity effect have focused on either the dissipative trans- 
port or the non-dissipative supercurrent. For low bi- 



ases redwithin the subgap transport region the dissipa- 
tive trans por t current is carried by multiple Andreev 
reflections [llj (MAR). MAR resonances occur when a 
sequence of Andreev reflections connect the high den- 
sity of states at the edge of the superconducting en- 
ergy gap. Sub-gap transport MAR resonances have been 
shown to be substantially altered by the sin gle electron 
states of a QD ||, EE E3 and the Kondo effect 0, El Ei | • 
Observation of the supercurrent in QD Josephson junc- 
tions presents a greater challenge which has been tack- 
led in a number of r ecent studies 0, EE El, El- In 
weakly coupled devices where the doublet state dom- 
inates the reversal of the supercurrent or 7r-junction 
has been demonstrated [H, 0, 17 1 in good agreement 
with theoretical predictions |d. ll8|T In the strongly cou- 
pled Kondo regime the supercurrent has been analyzed 
through current-biased V-I characteristics in QD Joseph- 
son i unctions [l6| and discussed in terms of the zero-bias 
peak in differential conductance measurements [19j with 
evidence of enhanced critical currents when tx > 1. In 
this report we study both the dissipative MAR trans- 
port and non-dissipative supercurrents, using a side-gate 
effect to electrically tune the scaling parameter tx- 



Devices were fabricated with a single uncapped InAs 
self- assembled QD with diameter and height of ~ 100 nm 
and ~ 20 nm, respectively. Conventional e-beam lithog- 
raphy and e-bcam evaporation techniques were used to 
deposit two Titanium/ Aluminium (5/100 nm) electrodes 
with a nanogap separation of less than 30 nm (Fig. Q] 
(a)). Additional devices fabrication detail can be found 
in references 0|23j|. We evaluate that A - 162 [ieV and 
T c ~ 1.1 K (Fig. Q] (a) and (d)). The back-gate is a de- 
generately Si doped GaAs layer buried 300 nm below the 
sample surface. The side-gate is placed about 200 nm 
away from the QD laterally. Transport measurements 



2 




a.) 



0.15 



Q.20 



0.25 



0.30 



0.35 



FIG. 1: (Color online)(a) Scanning electron microscope im- 
age of the sample studied and energy schematics for the re- 
order MAR (eVsd = A) and single quasi-particle tunneling 
resonances (eV s d = 2A). (b) Normal state stability diagram 
(B = 200 mT, V sg = V). Even (e) and odd (o) electron oc- 
cupations are indicated. White traces show the Kondo zero 
bias anomaly in the center of regions I and III. (c) Super- 
conducting state stability diagram (B — OmT, V sg = 0V). 
(d) Superconducting state G(V s d) traces for even electron oc- 
cupation Coulomb blockade (Vbg = 0.3 V, black trace), the 
center of region III ( Vbg — 0.335 V, red trace) and the center 
of region I (Vbg = 0.2 V, blue trace). 



were performed in a 3 He- 4 He dilution refrigerator with 
base temperature ~ 30 mK. The differential conductance 
was measured using conventional lock-in techniques with 
an ac excitation of V ac ~ 3 /iV. For measurement of the 
supcrcurrent a four-terminal setup with current-bias was 
applied. 

Fig. Q] (6) shows a false colorplot of the differen- 
tial conductance (G = dI/dV s d) taken by sweeping the 
back-gate Vb g and source-drain voltage V s d under an ap- 
plied magnetic field of 200 mT and zero side-gate bias 
(V sg = V). The applied field exceeds the Al lead criti- 
cal field (B c ~ 150 mT) such that the leads are in the 
normal state. In the Vb g range measured we identify 
three regimes with odd electron occupation, labeled as 
I, II, and III. The charging energy U is roughly esti- 
mated from the width of the Coulomb diamonds to be 
2.1,2.9 and 2.0 meV in regions I, II and III respectively. 
We evaluate the backgate leverarm, which relates Vb g to 
the energy in the QD0], to be a bg ~ 0.09 cV/V. Rc- 
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FIG. 2: (a) Normal state Coulomb peaks at V s d = 0V for a 
range of applied V sg . Curves from top to bottom are Vsg = 
+1.5 to —0.5 V in —0.5 V steps. Measurements are offset for 
clarity, (b) Summary of normal state (F) for regions I, II 
and III. (c) Summary (x) for regions I, II and III. (d) Tk 
evaluated in region III from the FWHM of the Kondo feature 
at B = 200 mT (o) and calculated from (r) and U ~ 2meV 
(•). The horizontal dashed line indicates A/fcs. The solid 
line indicates a linear best fit to the FWHM data. 



gions I and III exhibit Kondo zero-bias anomalies which 
are split and broadened by the Zecman energy respec- 
tively (g* = 6.1 ± 0.2 and 5.1 ± 0.3 for regions I and 
III respectively). In region III we estimate Tk from 
the full width at half maximum (FWHM= 2k B T K ) of 
Lorentzian fits, with a linear background subtraction, to 
the Kondo feature 2l| . By subtracting the Zeeman en- 
ergy from the FWHM a nominal Kondo temperature of 
Tk = 2.6 ± 0.3 K [tn ~ 1.4) is evaluated. In region I we 
observe the Kondo feature is already split at B = 200 mT 
indicating the Zeeman energy exceeds Tk , giving an up- 
per limit of T K < 0.82 ± 0.03 K (t K < 0.44). No Kondo 
feature is observed in region II indicating that Tk is much 
lower than the measurement temperature. 

When B = OmT, the leads are in the superconduct- 
ing state, Fig. Q] (c). In the even electron occupation 
regime where Coulomb blockade dominates (Fig. [H(rf)) 
prominent resonances at |el4<j| = 2 A are attributed to 
direct quasi-particle tunneling between the high density 
of states at the edge of the superconducting gap in the 
two leads. A weaker feature at \eV sc i\ = A is attributed 
to resonant single Andreev reflections or the l st -order 
MAR resonance (see Fig. [U (a)). In odd electron oc- 
cupation regions the spectrum of MAR features may be 
renormalized by the Kondo effect resulting in enhanced 



l st -order MAR features at the expense of single quasi- 



particle tunneling (12, 13, 14 1. Of the regions considered 
only region III displays enhancement of l s *-order MAR 
relative to the 2 A features, likely due to the higher Tk- 
A zero-bias conductance peak is observed in both even 
occupation regions and regions II and III, however this 
feature is absent in region I. In regions II and III we 
can eliminate the Kondo zero-bias anomaly as an origin 
of the superconducting state zero-bias peak because the 
magnetic field dependence does not show Zeeman split- 
ting. The zero-bias peak can therefore be re gar ded as a 
signature of supercurrent through the device [19(. A sig- 
nificant supercurrent feature is observed only in region 
III providing evidence of enhancement for high normal 
state Tr-Q. 

We will now focus on the effect of the side-gate on 
the normal state transport. Coulomb oscillations in the 
normal state for a range of V sg are shown in Fig. [5] 
(a). Coulomb peaks are shifted towards lower Vb g as 
V sg is increased with an evaluated leverarm of a sg ~ 
0.025af, 9 . We observe that V sg also alters the lineshape 
of the Coulomb peaks. Estimates of both the asym- 
metry in source (S) and drain (D) lead tunnel cou- 
plings (x = Ts.d/^d.s) and the total tunnel coupling 
(T = Ts + I\d) are obtained by fitting the even valley 
part of the Coulomb oscillation peaks using a Lorentzian 
expression following the method described in reference 
[ig | (see supporting information [20l| ) . Results of the av- 
erage tunnel coupling (T) and average tunnel coupling 
assymmetry (x) evaluated for the two Coulomb peaks in 
all three regions are plotted in Fig. \%l(b) and (c) respec- 
tively Values of (x) for regions I and III are found to 
be fairly constant while region II displays a decrease in 
(x) with increasing V sg which accounts for the increase in 
peak conductance of the fourth Coulomb peak observed 
in Fig. [2] (a). In all regions we observe no significant 
change in the g-factor or U when V sg is altered. We es- 
timate Tjc in experiment from the FWHM of the Kondo 
zero-bias anomaly[3| in region III for a range of V sg in 
Fig. [U (d). The dependence of Tk on V in the center of 
the odd electron occupation region is well formulated as 
T K ~ i\/!7rexp(-7r[//4r)[23]. The observed decrease 
in Tk calculated from the FWHM with increasing V sg is 
consistent with the decrease in that calculated from (r) 
evaluated from the Coulomb peaks in region III. These 
results indicate that in region III Tk may be tuned by 
controlling the tunnel coupling via the side-gate. In re- 
gions I and II (r) (and U) remains relatively constant 
when V sg is applied and we therefore observed no no- 
ticeable effect on Tk- T is determined by the effective 
overlap between lead states and the wavefunction of the 
confined electrons j23[. We expect that the side-gate mod- 
ulates the lateral confinement and displaces the electron 
wave function, resulting in a change in T. The influence 
of the side-gate is therefore highly dependent on specific 
symmetry of the orbital state and the corresponding ex- 
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FIG. 3: (Color online) Plots of G(V ad ,V bg ) displaying the 
subgap transport resonances in region III for V ag = —0.5 V 
(a) and V ag = +1.5 V (b). Horizontal lines indicate bias for 
quasi-particle tunneling resonance \eV 3d \ = 2A and l st -order 
MAR |eVsd| = A. Even (e) occupation regions are identified, 
(c) Plot of G(V s d) in the center of region III for a range of 
V S g. Traces from bottom to top are offset by +0.2e 2 /h for 
clarity. 



pansion of the wavefunction, which accounts for the dif- 
ferent behavior in regions I, II and III. We conclude that 
in region III the parameter Ik may be tuned in the range 
t K ~ 1.5 -> 0.8. 

We now focus on the superconducting transport in re- 
gion III as tK is tuned through unity. Details of the effect 
of V sg in regions I and II are included in the supplemental 
information [20j|. Fig. [21(a) shows a false color plot of dif- 
ferential conductance around region III for V sg = —0.5 V 
(a) and +1.5 V (6), focusing on the subgap transport res- 
onances described earlier. While the even occupation re- 
gions display similar subgap transport resonances when 
V sg = +1.5 and V sg = —0.5V, the odd occupation re- 
gion displays pronounced differences. The zero-bias con- 



ductance peak is observed for V SI 



-0.5 V but not for 



V sg = +1.5 V, indicating a suppression of the supercur- 
rent when V sg is increased. Fig. [3] (c) displays plots 
of G(V s d) at the center of the region III for a range of 
V S g. The magnitude of the zero-bias peak is gradually 
reduced and disappears at around V sg ~ +0.8 V. When 
the zero-bias peak diminishes a pair of small resonances 
are observed near V s( i = 0. These resonances may be at- 
tributed to higher order MAR. We also observe a distinct 
shift in the most prominent transport resonance from the 
l st -ordcr MAR feature at \eV s d\ = A to an intermediate 
feature between \eV s d\ = 2A and \eV s d\ = A. Note also 
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FIG. 4: (Color online) Plot of dV/dl for four terminal I- 
bias measurements in region III with V sg = —0.5 V (a) and 
V ag = +1.5 V (6). Dashed white lines indicate the trend of 
Isw- (c) Isw(Vb g ) in region III for a range of V sg . (d) I S w(ti<) 
for the center of region III. tx is estimated from the linear fit 
in Fig. Even (e) and odd (o) occupation is indicated. 



that the minima at \eV s d\ = 2 A are suppressed as V sg is 
increased. Viewed in the context of recent studies of the 
interplay between MAR and the Kondo effect II, 13. 14 1 
we attribute the pronounced change in the resonances to 
a reduction of the influence of Kondo effect on the spec- 
trum of MAR resonances. We note similarity between 
the subgap resonances in region III for high V sg (indicat- 
ing low Tk) with those observed in region I (shown in 
Fig. [I] (d)) where t K < 0.44. 

We analyze the transport further using a four ter- 
minal current bias measurement to study the super- 
current. We plot the four-terminal differential resis- 
tance (dV/dl) as a function of Vb g and current (/) for 
V, 



0.5 V and V sg = +1.5 V in Fig. g] (a) and (6) 
respectively. When measuring the non-dissipative cur- 
rent in the junction we must consider the dissipation in 
the surrounding circuit 24\ which determines the char- 
acteristics of the V-I curves. From consideration of the 
RCSJ model @, [H, 03, El we determine that in region 
III the junction is heavily overdamped [20I] with a non- 
zero resistance due to thermal phase diffusion [2^]. We 
qualitatively evaluate the relative evolution of the junc- 
tions intrinsic critical current using a 'switching current' 
(I S w), taken as the current value at the maximum dif- 
ferential resistance (white dashed lines in Fig. [2] (a and 
b)). Isw is influenced by the dissipation of the extrinsic 
circuit and is expected to be significantly reduced from 
the intrinsic critical current due to thermal fluctuations. 
In Fig. H](a) where V sg — —0.5 V, I sw is high in the odd 



occupation region and small but non-zero in the even 
occupation regions where Coulomb blockade dominates. 
Fig. [4] (c) shows the evolution of I sw (Vb g ) for a range 
of V sg . In good agreement with the V^-bias measure- 
ments we find that supercurrent feature is reduced to 
zero in the odd electron occupation region when V sg is 
high. Recent functional^ and numerical^, ,S] rcnormal- 
ization group studies predict a sharp drop in absolute I c 
when the device is tuned to the magnetic doublet regime. 
The overall trend of I sw (Vb g ) on V sg shown in Fig. [4] (d) 
therefore qualitatively matches the scenario of a phase 
transition between a Kondo dominated singlet state to 
a degenerate 'magnetic' doublet. In the magnetic dou- 
blet regime (tK < 1) the supercurrent is strongly sup- 
pressed by Coulomb blockade while in the Kondo singlet 
regime (tpc > 1) the supercurrent may be enhanced by 
the Kondo effect. The transition is observed at tx ~ IT 
in good agreement with the disappearance of the zero- 
bias peak in voltage bias measurements. Choi et aZ.[H 
used the numerical renormalization group to calculate I c 
as a function of tK and predicted the transition at tjc ~ 2. 
For tK > 2 a saturation of I c was predicted. Siano and 
Eggerjgj applied the Hirsch-Fye Monte Carlo method to 
predict a transition at tK ~ 1.1. Our experimental result 
is in good qualitative agreement with the predicted fea- 
tures in both of these studies. Discrepancies may arise 
from the non-ideal nature of real devices in which tunnel 
coupling is assymmetric and extrinsic environmental ef- 
fects may dominate. In some measurements [l7} the low 
supercurrent in the magnetic doublet (7r-junction) regime 
(when it « 1) has been measured, however in the de- 
vice considered here no supercurrent branch is observed 
implying that coherent processes are overcome by dissi- 
pation in the circuit and thermal fluctuations due to a 
poorly screened electromagnetic environment. 

In the present report the side-gate electrode is demon- 
strated to allow limited tuning of the device parameters 
in a region of constant electron number. This technique 
is used to elucidate the effect of Kondo correlations on 
both the supercurrent and subgap dissipative transport 
around the phase transition between magnetic double 
and Kondo singlet states. The onset of the transition 
was observed at tK — fesTW/A ~ 1.1. 
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